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摘 要 : 土地 利用 / 覆 被 变化 (Land Use and Cover Change , LUCC ) 是 全 球 环境 变化 的 重要 组 成 部 分 。 数 量 - 空 间 耦 
合 模型 因 能 有 效 模拟 土地 利用 斑 块 数量 及 其 在 空间 上 的 快速 变化 ,成 为 近年 来 土地 利用 研究 的 主要 手段 。 本 研究 


首先 对 比 斑 块 生成 士 地 利用 变化 模拟 模型 (Patch-Generating Land Use Simulation Model , PLUS) .未 来 土地 利用 模型 


(Future Land Use Simulation Model ,FLUS) 和 小 尺度 土地 利用 变化 及 其 空间 效益 模型 (Conversion of Land Use and Its 


Effects at Small Regional Extent,CLUE-S) 模 拟 2015 年 黑河 流域 中 游 甘 临 高 地 区 (张掖 市 甘 州 区 , 临 泽 县 和 高 台 县 ) 土 


地 利用 结构 ,筛选 出 最 适合 研究 区 的 空间 模拟 模型 ;然后 将 其 与 系统 动力 学 模型 (System Dynamics,SD ) 进 行 耦合 ， 
儿 成 数量 -空间 耦合 模型 ;最 后 利用 耦合 模型 预测 2030 年 研究 区 在 经 济 发 展 (Economic Development, ED) .生态 保 
护 (Ecological Protection ,EP) 和 协调 发 展 (Harmonious Development, HD ) 情 景 下 的 土地 利用 结构 ,并 进行 对 比 评价 。 
研究 结果 表明 :(1) 在 同一 像 元 尺度 下 ,PLUS 模型 的 质量 系数 (Figure of Merit,FoM) 高 于 FLUS 和 CLUE-S 模 型 ,上 且 3 


个 模型 的 数量 Kappa 系数 (K,,) 和 位 置 Kappa 系数 (Ko,) 从 高 至 低 为 :PLUS、CLUE-S 和 FLUS 模 型 ,说 明 PLUS 模 型 


在 本 研究 区 的 空间 拟 合 效 果 最 优 ; (2) PLUS-SD 耦合 模型 预测 的 不 同情 景 表 明 ,ED 情景 中 建设 用 地 和 耕地 迅速 扩 
张 但 森林 覆盖 率 较 低 ,EP 情景 中 的 森林 覆盖 率 明 显 提高 ,但 经 济 发 展 水 平 偏 低 ,HD 情景 中 城市 化 水 平 提 高 的 同时 


生态 环境 呈 良 性 发 展 。 本 研究 验证 了 PLUS-SD 耦合 模型 的 显著 优势 ,凸显 了 PLUS 模型 有 效 的 土地 利用 空间 模拟 


能 力 与 SD 模型 强大 的 数据 处 理 功 能 。 该 结果 可 为 土地 规划 政策 的 制定 提供 参考 。 


关键 词 : 黑河 流域 中 游 ， 耦 合 模型 ， 土 地 利用 / 奢 被 变化 ; 情景 预测 


20 世 纪 90 年 代 土地 利用 结构 变化 问题 开始 受 
到 关注 中 , 现 已 成 为 LUCC 研 究 的 热点 。LUCC 模 型 
是 研究 土地 利用 动态 变化 的 重要 手段 ,根据 各 种 模 
型 的 优势 可 以 分 为 数量 模拟 模型 .空间 模拟 模型 及 
耦合 模型 “”。 数 量 模拟 模型 主要 包括 SD (System 
Dynamics ) , GF(Gray Forecast) 和 MC (Markov Chain) 
等 模型 ,这 类 模型 能 有 效 地 模拟 和 预测 土地 利用 类 
型 的 数量 结构 ,但 难以 实现 空间 结构 模拟 “5 。 空 
间 模 拟 模 型 有 CA (Cellular Automaton) , CLUE- S 
(Conversion of Land Use and Its Effects at Small Re- 
gional Extent) , FLUS (Future Land Use Simulation 
Model) , ABM (Agent- Based Model) 和 PLUS ( Patch- 
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模式 。Liang 等 "9 基于 SD-CLUE-S 和 InVEST(Inter- 
grated Valuation of Environmental Services and Trad- 
eoffs ) 模 型 模拟 和 预测 了 2000—2018 年 张掖 绿洲 土 
地 利用 变化 对 像 元 扩 度 和 区 域 斥 度 碳 储 量 的 影 
响 。Ding 等 "建立 了 一 个 基于 FLUS 与 InVEST 模 
型 的 而 合 模型 ,有 效 地 模拟 了 东营 市 2030 年 不 同情 
景 的 人 居 环 境 质量 变化 。 

研究 表明 ,在 耦合 模型 中 ,空间 模拟 模型 的 性 
能 尤为 重要 ,而 不 同 模型 的 适用 区 域 存在 差异 。 
PLUS 模型 通过 随机 森林 算法 挖掘 地 类 转化 概率 , 善 
于 处 理 高 维 数据 "。CLUE-S 模 型 来 源 于 CLUE 模 
型 ,解决 了 了 CLUE 模型 不 适用 于 小 区 域 模拟 的 局 限 
性 * 中。FLUS 模 型 采用 神经 网 络 算法 获得 适宜 性 概 
率 , 有 利于 保证 模拟 精度 的 提升 ””。Peng 等 “将 
随机 森林 算法 与 CLUE-S 模 型 结合 ,充分 利用 随机 
森林 算法 在 处 理 高 维 共 线 性 问题 方面 的 优势 。 林 
丽 等 3 选择 LCM (Land Change Modeler) , FLUS CA- 
Markov 和 CLUFE-S 模 型 对 地 形 复杂 的 云南 山区 进行 
研究 ,最 终 发 现 不 同 土地 利用 类 型 的 适用 模型 存在 
差异 ,FLUS 模 型 更 适 于 模拟 建设 用 地 ,CLUE-S 模 型 
在 园地 存在 过 拟 合 现 象 。Liang 等 "对 武汉 市 进行 
研究 ,发 现 PLUS 模 型 在 该 区 域 的 模拟 性 能 优 于 
FLUS 模 型 。 由 此 可 见 , 不 同 空间 模拟 模型 在 不 同 区 
域 表 现 出 不 同 的 模拟 结 

位 于 中 国 西 北 干 量 区 的 张掖 市 甘 州 区 、 临 泽 县 
和 高 台 县 (本 文 简称 甘 临 高 ) ,处 于 河西 走廊 中 部 的 
黑河 流域 中 游 .青藏 高 原 和 蒙古 高 原 的 连接 处 , 具 
有 典型 的 区 域 代表 性 ,加 之 该 区 生态 环境 和 社会 经 
济 发 展 对 土地 利用 变化 较为 敏感 。 因 此 ,研究 该 区 
的 土地 利用 变化 具有 重要 意义 。 为 深入 探究 PLUS、 
CLUE-S 和 FLUS 空间 模拟 模型 在 该 区 的 适用 性 ,更 
好 地 预测 土地 利用 变化 ,笔者 选取 甘 临 高 地 区 为 研 
究 区 ,通过 对 比 PLUS FLUS 和 CLUE-S 模 型 模拟 结 
AR ,筛选 最 适 于 该 区 的 空间 模拟 模型 ;然后 将 其 与 
SD 数量 模拟 模型 耦合 , 预测 2030 年 不 同情 景 下 研 
究 区 的 土地 利用 格局 ,为 区 域 社会 经 济 可 持续 发 展 
提供 参考 。 


1 研究 区 概况 与 方法 


1.1 研究 区 概况 
黑河 是 中 国 第 二 大 内 陆 河 ,黑河 流域 中 游 的 张 


掖 市 甘 州 区 、 临 泽 县 和 高 台 县 的 是 本 文 的 研究 区 。 
该 区 地 理 位 置 为 98°57'~100°52’E、38°39’~39°59’'N 
(图 1)。 人 研究 区 的 气候 类 型 为 温带 大 陆 性 气候 ,年 
均 降 水 量 约 为 100~150 mm, 4E 3g 2€ AK Ek e à 
2047.9 mm ,年 均 气 温 在 8~9 % 区 间 波 动 。 研 究 区 地 
势 南北 高 .中间 低 ,主要 的 地 表 履 盖 类 型 为 耕地 E 
地 荒漠 等 的 。 该 区 位 于 干旱 区 ,生态 环境 脆弱 , 草 
地 的 萎缩 和 耕地 的 不 合理 扩张 加 剧 了 蕊 漠 化 。 
1.2 数据 来 源 及 处 理 

土地 利用 数据 、 自 然 环 境 数据 和 社会 经 济 数据 
均 来 源 于 中 国 科 学 院 资源 环境 数据 共享 中 心 (http:// 
www.resdc.cn)。 除 气温 、 降 水 、 人 口 和 GDP 数 据 的 
空间 分 辨 率 为 1 km, 其余 数据 的 空间 分 辩 率 为 30 m 
( 表 1)。 本 研究 将 所 有 数据 的 空间 分 辨 率 重 采样 为 
30 m。 不 同 像 元 与 道路 河流、 居民 点 的 距离 数据 
均 是 ArcGIS 10.2 软 件 中 基于 相应 矢量 数据 计算 欧 
REKHA 
13 模型 与 方法 
1.3.1 土地 利用 空间 格局 模拟 及 验证 

(1) 3 种 空间 模型 介绍 

PLUS 模型 包括 基于 土地 扩张 分 析 策 略 的 转化 
规则 挖掘 框架 和 基于 多 类 型 随机 斑 块 种 子 机 制 的 
CA 模型 两 大 模块 (图 2)。CLUE-S 模 型 包括 非 空 
间 土 地 需求 模块 和 空间 分 配 模 块 两 大 部 分 *”。 本 
人 研究 采用 线性 插值 法 计算 2000 一 2015 年 间 不 同 地 
类 的 需求 面积 ,不 同 地 类 与 驱动 因素 之 间 的 关系 数 
据 来 源 于 Logistics 回归 方程 。 研 究 采 用 ROC (Rela- 
tive Operating Characteristics ) 方法 检验 Logistics 回 
归 的 结果 , 当 ROC>0.7, 说 明 回 归 方 程 能 够 较 好 地 解 
释 地 类 结构 5, 然后 采用 空间 分 配 模 块 进行 空间 结 
构 模 拟 。FLUS 模型 由 2 个 模块 构成 ”” ,第 一 个 模块 
采用 神经 网 络 算法 基于 初始 年 份 的 土地 利用 数据 
与 驱动 因素 栅 格 数据 获得 适宜 性 概率 ,第 二 个 模块 
基于 轮 盘 赌 法 则 的 自 适 应 惯性 苋 争 机 制 研究 不 同 
驱动 因素 影响 下 各 地 类 的 转化 。 

(2) 模型 精度 验证 指标 

主要 采用 3 种 模型 精度 验证 方法 。 第 一 种 方法 
是 逐 像 元 精度 评价 方法 ,评价 指标 为 FoM2 2 。 第 二 
种 方法 是 空间 格局 比较 法 ,评价 指标 为 景观 格局 相 
似 度 ,本 研究 采用 Fragstats 4.2.1 软件 计算 模拟 与 实 
际 的 土地 利用 结构 图 的 景观 指数 ,进而 获得 景观 格 
局 相似 度 Sim" ,公式 如 下 : 
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图 采用 自然 资源 部 标准 地 图 制作 , 审 图 号 GS(2019)3333 号 ,对 底 图 边界 无 修改 。 下 同 。 


1 研究 区 在 甘肃 省 的 地 理 位 置 (a) 、 高 程 (b)、 地 貌 类 型 (c) 和 2015 年 土地 利用 空间 分 布 (d) 


Fig. 1 Geographical location of the study area in Gansu Province (a), elevation (b), landform type (c) and the spatial distribution 


map of land use in 2015 (d) 


表 1 研究 数据 源 信息 


Tab.1 Data source information of the research 


土地 利用 数据 土地 利用 数据 (2000 年 ,2015 年) 30m http://www.resde.cn 
然 环境 数据 年 均 气 温 (2015 年 ) 1 km 
年 均 降 水 (2015 年 ) 1 km 
高 程 30m 
坡度 30m 
3i qn) 30m 
社会 经 济 数据 人 口 (2015 年 ) 1 km 
GDP(2015 4E.) 1 km 
像 元 与 道路 的 距离 数据 30m 
像 元 与 河流 的 距离 数据 30 m 
像 元 与 居民 点 的 距离 数据 30m 


"NEL. 
so 


L-L 
| CX 2 aus (1) 
SUP in Sean ti EVA /表示 模拟 土地 利用 结构 图 
的 景观 指数 ; L, 表示 实际 土地 利用 结构 图 的 景观 指 
数 。 第 三 种 方法 是 Kappa 系数 ,该 系数 主要 被 用 于 
验证 空间 模拟 精度 。 本 文 使 用 PONTIUS 发 展 的 


Kappa 系 数量 化 位 置 错误 和 数量 错误 呈 ,采用 IDRI- 
SI 17.0 软件 计算 数量 Kappa 系数 (K,,) f Er Kappa 
系数 (Ks) 和 标准 Kappa 系数 (Kw ) o 
1.3.2 基于 SD 模型 的 未 来 情景 预测 

(1) SD 模型 介绍 

SD 模型 于 1956 年 建立 ,该 模型 是 生态 环境 和 
经 济 发 展 复合 系统 的 仿真 模拟 实验 室 ,在 数据 模拟 
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图 2 本 研究 流程 


Fig.2 Flow chart of this study 


方面 功能 强大 "™“。 本 人 研究 使 用 的 SD 模型 包括 人 
口 .经济 .生产 力 和 气候 4 个 子 系统 (图 3)。 甘 临 高 
地 区 的 各 变量 数据 参考 张掖 地 区 ,各 系统 变量 的 数 
值 来 源 于 《甘肃 发 展 年 鉴 兴 张掖 统计 年 鉴 》 和 《中 
城市 统计 年 鉴 } 等 文献 。 使 用 SPSS 23 软件 进行 
统计 分 析 ,获得 不 同 变量 间 的 数量 关系 ,SD 模型 模 
拟 的 时 间 范 围 为 2000 一 2030 年 ,时 间 步 长 为 1 a。 

(2) SD 模型 检验 

SD 模型 历史 仿真 检验 的 相对 误差 计算 公式 为 : 
a | 模拟 值 -历史 值 | 
相对 误差 = Wü 

为 验证 模型 仿真 效果 ,本文 确定 历史 检验 时 间 
为 2015 年 ,历史 检验 的 相对 误差 均 低 于 5$% ,这 说 明 
模型 的 模拟 精度 较 高 ,能 够 正确 预测 土地 系统 结构 
( 表 2)。 

(3) 情景 选择 

“十 二 五 ?规划 期 间 张 掖 市 CDP 年 均 增 速 为 
10.596 ,人 口 年 均 增 长 率 为 5.2%o， 十 三 五 "规划 期 
间 GDP 年 均 增 长 率 为 8.5% ,综合 考虑 将 HD ED 和 
EP 情景 下 的 GDP 增 长 率 分 别 设置 为 12%、14% 和 


x100% (2) 


9%( 表 3); 由 于 二 孩 政策 的 实施 ,本 人 研究 将 HD ED 
和 EP 情景 下 的 人 口 增长 率 分 别 设置 为 6%o、10%o 和 
4%o。 甘 肃 省 城镇 体系 规划 (2013 一 2030 年 ) 中 要 求 
规划 期 末 张 掖 绿洲 生态 城乡 统筹 发 展区 城镇 化 率 
应 达到 61.8%, 据 此 本 研究 将 HD、ED 和 EP 情景 下 
的 城镇 化 率 分 别 设置 为 61.8% .65% 和 60%。IPCC 
报告 指出 2016 一 2035 年 全 球 气温 和 降水 会 逐年 增 
In, Bea AES” ,HD ED 和 EP 情景 下 的 
年 均 气 温 与 降水 变化 量 都 设置 为 正 值 。 根 据 统计 
年 鉴 中 的 技术 研发 投入 数据 确定 技术 参数 。 


2 结果 与 分 析 


2.1 PLUS 模型 与 其 他 模型 的 模拟 结果 

为 对 比分 析 CLUE-S .PLUS 和 FLUS 模 型 的 模拟 
效果 ,通过 重 采样 将 研究 像 元 尺度 统一 为 50 m。 当 
模拟 尺度 为 50 m 时 ,耕地 林地、 农村 建设 用 地 和 未 
利用 地 的 ROC 最 高 ( 表 4)。 由 此 可 知 ,CLUE-S 模 型 
在 本 研究 区 的 最 佳 研究 尺度 为 30 m。PLUS 模 型 在 
30 m,50 m,100 m 和 200 m HHK I RE F KY FoM 分别 
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< 耕地 > 


<Time> 


第 一 产业 固定 
资产 投资 系数 


EKES 


降水 后 一 一 一 一 <Time> a 


< 第 一 产业 
司 定 资产 
Hb e 


~ 
UT 
畜牧 业 投 资 系 数 


水 产 投资 系数 


畜牧 业 投资 


温 > <Time> 


一 产业 固定 T | 
农业 投资 系数 资产 投资 额 > 技术 进步 
图 3 SD 模型 结构 框架 
Fig.3 Flow chart of SD model 
表 2 SD 模型 模拟 结果 的 历史 性 检验 
Tab.2 Historical test of SD model' s simulation results /hm 
检验 耕地 林地 草地 水 域 城市 建设 用 地 农村 建设 用 地 未 利用 地 

历史 值 (2015 年 ) 217147 11625 149341 24752 6974 15450 642389 

模拟 值 (2015 年 ) 215784 11554 149322 24695 6938 14736 644650 

相对 误差 /% 0.63 0.61 0.01 0.23 0.51 4.62 0.35 

表 3 不 同情 景 的 参数 设 定 
Tab.3 Parameter setting of different scenarios 

变量 人 口 增长 率 /‰o GDP 增 长 率 /% 年 降雨 变化 /(mm*a') 年 均 温 变化 /% .ar ) ”技术 系数 城镇 化 率 
HD 情景 6 12 0.350 0.020 400 43% 线 性 增 至 61.8% 
ED 情景 10 14 0.550 0.030 500 43% 线 性 增 至 65% 
EP 情景 4 9 0.15 0.010 500 43% 线 性 增 至 60% 

表 4 CLUE-S 模 型 的 ROC 指 数 
Tab.4 ROC index of CLUE-S model 

分 辨 率 /m 耕地 林地 草地 水 域 城市 建设 用 地 农村 建设 用 地 未 利用 地 

30 0.731 0.566 0.57 0.509 0.569 0.684 0.758 

50 0.935 0.696 0.737 0.722 0.915 0.834 0.893 

100 0.935 0.698 0.729 0.728 0.867 0.832 0.892 
200 0.935 0.68 0.738 0.72 0.928 0.833 0.892 
为 0.234、0.184、0.145 和 0.152, 相同 像 元 尺度 下 ”上 度 均 为 30 m。 


FLUS 模 型 的 模拟 精度 分 别 为 0.147、0.133、0.128 和 
0.117。 据 此 可 得 ,PLUS 和 FLUS 模 型 的 最 佳 拟 合 尺 


在 50 m WRR E F PLUS, FLUS 和 CLUE-S 模 
型 的 FoM 值 分 别 为 0.184、0.133 和 0.133 ,PLUS 模型 
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在 模拟 精度 上 优 于 FLUS 和 CLUE-S 模 型 。PLUS、 
CLUE-S 和 FLUS 模 型 的 K, 分 别 为 0.9625 .0.9444 和 
0.9188 , Ku. 分 别 为 0.9456、0.9193 Fil 0.8818 , Kis 
分 别 为 0.9455 .0.9191 和 0.8818。3 个 模型 的 Kappa 
系数 值 均 大 于 0.8, 且 玉 , 高 于 Km 说明 3 个 模型 的 
数量 预测 能 力 均 优 于 位 置 预测 能 力 。PLUS 模 型 的 
3 种 Kappa 系 数 最 高 , 且 均 大 于 0.95 ,说 明 模 型 模拟 
结果 与 实际 情况 基本 一 致 。 

本 文采 用 15 个 景观 指数 比较 3 个 模型 在 挖 所 
景观 变化 方面 的 能 力 。Fragstats 4.2.1 软件 计算 的 结 
果 表 明 ( 表 5),PLUS 模 型 有 9 个 景观 指数 更 贴近 实 
际 的 景观 指数 :NP、PARA_MN PARA_MD, PA- 
RA, RA, PARA. CV, ENN MN, ENN RA, ENN MD 
和 ENN_SD;CLUE-S 模 型 有 4 个 景观 指数 与 实际 的 
景观 指数 相符 : PARA, AM, PLADJ, ENN. MD 和 
ENN. CV;FLUS #4! A LPI, PARA, SD 和 ENN_AM 3 
个 景观 指数 值 与 实际 的 景观 指数 相近 。PLUS、 
FLUS 和 CLUE-S 模 型 的 景观 相似 度 Sim 分 别 为 
0.64、0.21 和 0.46, 这 说 明 PLUS 模 型 在 景观 格局 模 
拟 方面 优势 明显 。 

总 体 上 ,PLUS 模型 的 整体 模拟 效果 明显 优 于 


CLUE-S 和 FLUS 模 型 ,CLUE-S 模 型 在 地 势 平坦 地 区 
模拟 准确 度 高 于 FLUS 模 型 ,而 FLUS 模 型 在 地 势 起 
伏 较 大 地 区 的 模拟 效果 优 于 CLUE-S 模 型 (图 4)。 
(1) 在 临 泽 县 、 高 台 县 和 甘 州 区 ,CLUE-S 模 型 模拟 
的 地 形 平坦 处 的 耕地 范围 略 大 于 实际 范围 ,但 是 准 
确 度 较 高 。 在 高 台 县 南部 海拔 高 处 ,CLUE-S 模 型 模 
拟 的 草地 范围 大 于 实际 范围 ;在 临 泽 县 和 甘 州 区 的 
北部 海拔 高 处 ,CLUE-S 模 型 模拟 的 草地 范围 小 于 实 
际 范 围 ; 在 地 形 平坦 的 中 部 地 区 , CLUE-S 模型 模拟 
的 建设 用 地 范围 略 小 于 实际 范围 ,在 海拔 高 处 ,该 
模型 模拟 的 林地 和 建设 用 地 范围 大 于 实际 范围 。 
(2) 在 临 泽 具 和 甘 州 区 ,FLUS 模型 模拟 的 耕地 范围 
显著 大 于 实际 范围 ,在 高 台 县 FLUS 模型 模拟 的 耕 
地 范围 小 于 实际 范围 ;FLUS 模 型 模拟 的 海拔 高 处 的 
草地 范围 与 实际 情况 较为 接近 ,模拟 精度 较 CLUE-S 
模型 高 ,但 略 低 于 PLUS 模型 ;FLUS 模 型 模拟 的 林地 
和 建设 用 地 范 团 较 实 际 范围 略 大 。(3) PLUS 模型 的 
模拟 准确 度 显 著 高 于 CLUE-S 和 FLUS 模 型 ,其 模拟 
结果 基本 与 实际 范围 相符 。 与 FLUS 模型 相 比 ， 
CLUE-S 和 PLUS 模型 模拟 的 农村 建设 用 地 范围 与 实 
际 更 相符 。 甘 州 区 , 临 泽 县 和 高 台 县 的 城市 建设 用 


表 5 不 同 模型 2015 年 模拟 结果 的 景观 指数 


Tab.5 Landscape data of simulation results of different models in 2015 


景观 指数 PLUS(50 m) FLUS(50 m) CLUE-S(50 m) 实际 (50 m) 
NP (Number of patches) 7122.00a 20303.00c 11957.00b 4182.00 
LPI (Largest patch index )/% 43.48c 52.71a 53.06b 50.89 
PARA_MN (Mean perimeter-area ratio ) 456.47a 613.45c 536.94b 231.18 
PARA. AM (Area-weighted mean perimeter-area ratio ) 36.94b 63.30c 36.17a 33.58 
PARA. MD (Median perimeter-area ratio ) 366.67a 640.00c 600.00b 22222 
PARA. RA (Perimeter-area ratio range ) 7790.01a 716.04c 790.70b 789.71 
PARA. SD (Perimeter-area ratio standard deviation ) 266.51c 206.07a 250.83b 123.03 
PARA. CV (Perimeter-area ratio coefficient of variation ) 58.38a 33.59c 46.72b 53.22 
ENN MN (Mean euclidean nearest neighbor distance )/m 276.76a 150.13c 233.80b 464.08 
ENN AM (Area- weighted mean euclidean nearest neigh- 183.61b 146.66a 131.66c 159.72 
bor distance )/m 
ENN. MD (Median euclidean nearest neighbor distance )/m 141.42a 111.80c 141.42a 223.61 
ENN RA (Euclidean nearest neighbor distance range )/m 26382.45a 20753.12c 24047.88b 34894.36 
ENN SD (Euclidean nearest neighbor distance standard 751.54a 320.79c 516.44b 1087.47 
deviation )/m 
ENN. CV (Euclidean nearest neighbor distance coefficient 271.55c 213.68b 220.89a 234.33 
of variation )/m 
PLADJ (Percentage of like adjacencies )/% 95.38b 92.09c 95.48a 95.80 
Sim 0.64 0.21 0.46 


注 :a 指 代 模拟 图 像 与 实际 图 像 的 景观 指数 最 相符 ;b 次 之 ;c 表 示 最 不 相符 。 
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图 4 不 同 模型 2015 年 的 模拟 结果 对 比 


Fig. 4 Comparison of simulation results of different models in 2015 
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地 发 展 过 程 中 均 出 现 了 “ 飞 地 式 ” 扩 张 现象 ,3 个 模 
型 在 模拟 这 类 建设 用 地 时 均 失 效 , 这 或 许 是 驱动 因 
素 与 地 类 变化 之 间 时 间 滞 后 的 结果 。 
2.2 PLUS-SD 耦合 模型 的 不 同情 景 预测 结果 

HN EC , 临 泽 县 和 高 台 县 各 地 类 变化 的 区 域 多 
位 于 地 类 边缘 ,以 边缘 式 扩张 为 主 (图 $, 表 6~ 表 
7)。(1) 基于 30m 分 辩 率 土地 利用 数据 模拟 的 结果 
表明 ,EP 和 HD 人 情景 中 2015 一 2030 年 期 间 耕 地 面积 
增加 ,其 来 源 均 有 草地 。EP 情 景 下 耕地 面积 减少 ， 
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主要 转化 为 林地 和 草地 ,说 明 政 府 采 取 退 耕 还 林 还 
草 措 施 保护 生态 环境 。(2) ED 情景 中 城市 化 水 平 显 
著 提高 ,技术 发 展 迅 速 ,2030 年 林地 面积 逐渐 减少 ， 
主要 转化 为 耕地 。EP 和 HD 情景 林地 面积 增加 ,HD 
情景 中 部 分 未 利用 地 转化 为 林地 ,EP 情景 增加 的 林 
地 主要 来 自 耕 地 。(3) 2030 年 ED 和 HD 情景 的 草地 
面积 低 于 2015 年 ,EP 情景 则 相反 。(4) 城市 发 展 导 
致 热岛 效应 和 雨 岛 效 应 ,全球 变 暧 促使 降水 增加 ， 
模拟 结果 显示 3 种 情景 的 水 域 面积 均 增 加 。(5) 


99°0'E 99?20'E 99?40'E 100?0'E 


z 8 Zz 
$ £7 (CD HD 情 景 -高 台 E 
a Y A + 
和 " 
4 Zz 
© [=] 
= = 
z ES 
© 
全 E 
S$ £ £ 
$ S 
z£ £ 
$5 3 
99°0'E 99°20'E 99?40'E 100?0'E 
100?40'E 9990'E 99°20'E 99?40'E 100°0'E 
E E c 
+ & T (c2) EP 情 景 -高 台 县 Z 
a x S 
e ^ 
Z Zz 
© © 
= = 
z 8 $ 
© 
-- z 
$ £ E 
$ $ 
z Zz 
£ $ S 
a a a 
99°0'E 99?20'E 99°40'E 100°0'E 
100°40'E 99°0’E 99?20'E 99?40'E 100°0'E 
Z £ Z 
* & | (3) ED 情景 -高 台 县 Z 
人 "t + 
e ss 
区 区 
© Q 
= = 
z 8 $ 
[-] 
NE- z 
$ s x 
$ EA 
z § P 
$5 $ 


99°0'E 99?20'E 99?40'E 100?0'E 


图 例 ma 耕地 mmm 林地 ES 草地 Eke 一 | 城市 建设 用 地 L3 农村 建设 用 地 号 未 利用 地 


图 5 不 同情 景 中 的 2030 年 土地 利用 结构 


Fig. 5 Land use structure in different scenarios in 2030 
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表 6 不 同情 景 的 未 来 各 土地 利用 类 型 面积 预测 结果 


Tab.6 Future prediction results of land use types under different scenarios /hm 
"T er SD 模型 (2030 4E) 

ED 情景 EP 情景 HD 情景 
耕地 217147 261740 203469 225155 
林地 11625 7958 14747 12165 
草地 149341 137645 152568 146553 
水 域 24752 28097 25413 26704 
城市 建设 用 地 6974 14253 12070 13166 
农村 建设 用 地 15450 14732 14845 14690 
未 利用 地 642389 603253 644568 629247 


表 7 不 同情 景 的 未 来 各 土地 利用 转移 和 矩阵 


Tab.7 Future land use transition matrix under different scenarios /hm 
情景 2015 年 /2030 年 耕地 林地 草地 水 域 ”城市 建设 用 地 农村 建设 用 地 未 利用 地 总 计 
ED 情景 ”耕地 215340 18 272 239 292 0 987 217147 
林地 2803 7648 421 112 68 0 573 11625 
草地 6215 82 134401 435 283 0 7924 149341 
水 域 2825 27 718 17973 150 0 3059 24752 
城市 建设 用 地 414 7 16 5 6430 0 103 6974 
农村 建设 用 地 708 0 0 9 14732 0 15450 
未 利用 地 33435 175 1818 9332 7022 0 590606 642389 
总 计 261740 7958 137645 28097 14253 14732 603253 1067678 
EP 情景 ”耕地 184921 6293 11703 2506 1943 26 9754 217147 
林地 1341 6718 2020 150 106 0 1289 11625 
草地 4176 730 121801 1467 747 0 20419 149341 
水 域 1271 112 1856 18428 99 1 2985 24752 
城市 建设 用 地 246 32 94 5 6354 1 242 6974 
农村 建设 用 地 491 9 19 5 50 14813 64 15450 
未 利用 地 11022 854 15074 2851 2770 3 609815 642389 
总 计 203469 14747 152568 25413 12070 14845 644568 1067678 
HD 人 情景 ” 耕地 217015 2 4 17 25 0 83 217147 
林地 139 10715 11 133 54 0 572 11625 
草地 2181 48 146489 188 94 0 341 149341 
水 域 176 60 17 21684 118 0 2697 24752 
城市 建设 用 地 22 6 1 1 6874 0 70 6974 
农村 建设 用 地 748 1 0 1 12 14690 0 15450 
未 利用 地 4874 1334 31 4679 5989 0 625483 642389 
总 计 225155 12165 146553 26704 13166 14690 629247 1067678 


2030 年 3 种 情景 的 城市 建设 用 地 面积 均 增 加 ,主要 "- 
来 自 于 未 利用 地 。ED .EP 和 HD 情景 的 农村 建设 用 3 讨论 
地 均 呈 不 同 程度 下 降 , 可 能 是 大 量 农村 人 口 向 城 
人 口 转变 ,(6) 3 种 情景 的 未 利用 地 主要 向 耕地 和 城 31 PLUS-SD 模 型 在 甘 临 高 地 区 的 适用 性 分 析 
市 建设 用 地 转化 ,ED 情景 转化 面积 最 大 ,HD 情景 次 (1) 通过 PLUS FLUS 和 CLUE-S 模 型 在 甘 临 高 
之 ,EP 情景 转化 面积 最 小 。 地 区 的 模拟 结果 ,可 以 发 现 PLUS 模 型 在 位 置 预 测 、 


oy 
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数量 预测 和 景观 格局 相似 度 方面 均 优 于 FLUS 和 
CLUE-S 模 型 ,CLUE-S 模 型 在 地 形 较 为 平坦 地 区 模 
拟 精 度 较 高 ,而 FLUS 模型 在 海拔 高 处 模拟 效果 略 
优 于 CLUE-S 模 型 。 前 人 的 研究 得 出 了 类 似 结论 ， 
Liang 等 ”的 研究 表明 ,PLUS 模 型 在 模拟 历史 土地 
利用 变化 过 程 方面 优 于 FLUS 模 型 ;济南 市 区 土地 
利用 变化 模拟 结果 显示 FLUS 模 型 模拟 结果 的 Kap- 
pa 2 AUKAR F CLUE-S Bi AY 07 , (2) 本 研究 表明 ， 
FLUS 45 PLUS 模型 的 研究 格 网 愈 密集 ,模拟 结果 愈 
优 ,而 CLUE-S 模 型 并 未 随 着 格 网 密度 增加 而 得 到 
更 优 的 模拟 结果 。 张 永 民 等 3 基于 Logistic 回归 模 
型 采用 100 m、250 m、500 m,750 m Fil 1000 m 这 5 种 
栅 格 尺度 研究 土地 利用 与 不 同 影响 因子 的 关系 , 结 
果 显 示 500 m 尺 度 下 二 者 的 相关 性 最 高 。CLUE-S 
模型 基于 传统 的 二 值 逻辑 斯 蒂 回 归 确 定 地 类 空间 
结构 与 影响 因素 之 间 的 相关 关系 ,其 研究 尺度 不 宜 
过 度 密 集 ; 而 PLUS 模型 在 模拟 细 尺 度 元 胞 变化 方 
面 能 力 较 强 ;FLUS 模型 通过 增强 空间 分 配 和 需求 模 
拟 之 间 的 相互 映射 以 达到 ”* 紧 耦合 的 目标 "。(3) SD 
模型 适用 于 模拟 土地 结构 中 不 同 影响 因素 之 间 的 
数量 关系 ,通过 情景 分 析 模拟 政策 变化 的 效果 。 
本 研究 将 SD 与 PLUS 模型 耦合 ,获得 了 不 同情 景 下 
的 土地 利用 结构 ,有 利于 生态 和 经 济 的 可 持续 发 
展 。 不 少 研究 对 未 来 不 同情 景 下 的 土地 利用 格局 
YET. Lu SE "3€ SD fIICCDM (Coupling Co- 
ordination Degree Model) Z4 ,对 武汉 市 经 济 -资源 - 
环境 系统 的 模拟 结果 进行 了 分 析 和 评价 ,清晰 地 揭 
示 了 经 济 -资源 -环境 系统 中 各 因素 的 耦合 关系 。 
Wu 55$ y SD-CLUE-S HATEN ,发 现 综合 性 模 
型 能 够 更 好 地 模拟 上 海宝 山地 区 景观 生态 系统 价 
值 的 动态 变化 。 
3.2 不 同情 景 模 拟 结果 的 对 比分 析 

甘 临 高 地 区 的 地 表 环 境 较 恶 劣 ,沙漠 和 戈壁 面 
积 较 大 。 本 人 研究 的 ED 情景 和 EP 情景 偏重 于 经 济 发 
展 或 环境 保护 ,不 利于 生态 和 经 济 协调 发 展 ;HD f 
景 兼 具 生态 友好 和 经 济 建设 的 优点 。 前 人 研究 表 
8j] , 甘 临 高 地 区 存在 草地 荒漠 化 现象 ,草地 生态 环 
境 脆 弱 ,耕地 对 水 分 和 养分 的 需求 量 大 ,不 合理 地 
开垦 耕地 将 加 剧 土地 荒漠 化 进程 “。2030 年 的 HD 
情景 中 有 少量 草地 和 农村 建设 用 地 转化 为 耕地 , 耕 
地 主要 来 自 于 未 利用 地 ;退耕 还 林 的 加 大 ,对 未 利 
用 地 的 开发 力度 促使 林地 面积 增加 ,有 利于 发 挥 国 


土 空间 的 生态 功能 。HD 情景 中 城市 扩张 的 主要 方 
向 也 是 未 利用 地 , 随 着 技术 的 发 展 ,未 利用 地 的 利用 
率 和 转化 率 显著 提升 。 近 年 来 高 台 等 地 的 绿洲 扩张 
明显 , 且 绿 洲 波动 带 多 位 于 绿洲 和 荒漠 交界 处 5 。 
随 着 城镇 化 进程 的 加 快 会 导致 建设 用 地 扩展 ,2033 
年 建设 用 地 面积 将 高 达 4948.90x10: hm? , TE HD 
情景 中 ,政府 需 严 格 管控 建设 用 地 ,防止 建设 用 地 过 
度 占 用 耕地 ,保护 永久 性 基本 农田 ,要 集约 且 合 理 地 
安排 建设 用 地 ,形成 科学 有 序 的 国土 空间 开发 格 
E, 


4 结论 


vd 


本 研究 发 现 PLUS 模型 在 景观 格局 .位 置 和 数 
量 模拟 方面 的 精度 均 优 于 CLUE-S 和 FLUS 模型 , 
CLUE-S 模 型 在 研究 区 地 势 低 平地 区 模拟 精度 较 高 ， 
而 FLUS 模 型 在 海拔 高 处 模拟 效果 较 CLUE-S 优 ;对 
比分 析 PLUS-SD 耦合 模型 模拟 的 甘 临 高 地 区 2030 
年 ED、EP 和 HD 情景 下 的 土地 利用 情况 ,发 现 3 种 
情景 在 土地 利用 类 型 变化 速率 方面 存在 差异 ;ED 情 
景 中 土地 利用 类 型 变化 速率 最 快 ,经 济 迅猛 发 展 ， 
但 是 环境 破坏 严重 ,林地 大 面积 减少 ,EP 情景 中 土 
地 利用 类 型 变化 速率 偏 慢 , 存 在 明显 的 退耕 还 林 还 
WMA (AA AR , HD 情景 中 建设 用 地 占用 
耕地 的 现象 较 少 ,经 济 平 稳 运 行 ,有 利于 社会 可 持 
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Land use in the Gan-Lin-Gao region of middle reaches of Heihe River Basin 
based on a PLUS-SD coupling model 


JIANG Xiaofang", DUAN Hanchen'*, LIAO Jie", SONG Xiang", XUE Xian™’ 
(1. Key Laboratory of Desert and Desertification, Northwest Institute of Eco-Environment and Resources, Chinese 
Academy of Sciences, Lanzhou 730000, Gansu, China; 2. University of Chinese Academy of Sciences, Beijing 
100049, China; 3. Drylands Salinization Research Station, Northwest Institute of Eco-Environment and Resources, 
Chinese Academy of Sciences, Lanzhou 730000, Gansu, China) 


Abstract: A quantity-space coupling model can effectively simulate the number of land-use patches and their rapid 
spatial changes. Thus, this model has become the main method of land use research in recent years. This research 
first compared the simulated land-use structure of the patch-generating land-use simulation model (PLUS), future 
land-use simulation model (FLUS), and conversion of land use and its effects at small regional extent (CLUE-S) 
in the Gan- Lin-Gao region (Ganzhou District, Linze County and Gaotai County of Zhangye City) of the middle 
reaches of the Heihe River Basin in 2015. In addition, this research selected the suitable spatial simulation model 
for the study area; then, it was coupled with the system dynamics model (SD) to form a quantitative-spatial coupling 
model. Finally, the coupling model was used to predict the land-use structure of the economic development (ED), 
ecological protection (EP), and harmonious development (HD) scenarios, and comparative evaluation was 
conducted. The research results showed that (1) at the same pixel scale, the figure of merit of the PLUS model 
was higher than that of the FLUS and CLUE-S models, and the K,, and Ki, of the three models were PLUS, 
CLUE-S, and FLUS from high to low, indicating that the PLUS model had the best spatial fitting effect in this 
study area. (2) Moreover, the different scenarios predicted by the PLUS- SD coupling model showed that the 
construction land and arable land in the ED scenario expanded rapidly, but the forest coverage rate was low. 
Furthermore, the forest coverage rate in the EP scenario increased significantly, but the level of ED was low, and 
the urbanization in the HD scenario improved when the ecological environment was gradually developing. This 
research determined the significant advantages of the PLUS- SD coupling model and highlighted the effective 
land-use spatial simulation capabilities of the PLUS model and powerful data processing capabilities of the SD 
model. The results can provide reference for the formulation of land planning policies. 
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